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1^ ' Abstract. We calculate the phase diagrams for a harmonically trapped ideal 

gas mixture of fermionic atoms and bosonic molecules in chemical and thermal 
equilibrium, where the internal energy of the molecules can be adjusted relative 
to that of the atoms by use of a tunable Feshbach resonance. We plot the 

C^ ■ molecule fraction and the fraction of Bose-condensed molecules as functions of 

ryp^ ' the temperature and internal molecular energy. We show the paths traversed 

in the phase diagrams when the molecular energy is varied either suddenly 

jS I or adiabatically. Our model calculation helps to interpret the adiabatic phase 

e, diagrams obtained in recent experiments on the BEC-BCS crossover, in which 

. the condensate fraction is plotted as a function of the initial temperature of the 

I ' Fermi gas measured before a sweep of the magnetic field through the resonance 

y ' region. 
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Figure 1. Illustration of the atom-molecule phase diagram for a Fermi gas with 
resonant interactions. Starting at some initial point in the phase diagram such as 
A, the system can be continuously transformed from a gas of fermionic atoms to a 
gas of bosonic molecules at some point B by adjusting the resonant state energy. 



1. Introduction 

In recent experiments, a Feshbacli resonance is utilized to explore the crossover from 
Bardeen-Cooper-SchriefFer (BCS) superfluidity in a dilute Fermi gas to Bose-Einstein 
condensation (BEC) in a dilute gas of diatomic molecules PI^IISIEI- A qualitative 
sketch of the phase diagram for this system is shown in Figure ^ By tuning the 
resonant state energy, the system can be continuously transformed from a gas of 
fermionic atoms, starting for example at point A and moving along some path to a 
gas of bosonic diatomic molecules at point B. The specific path taken depends on the 
rate at which the sweep is performed. 

In this article, we study the two limiting cases of an adiabatic sweep or a sudden 
change in order to gain a qualitative understanding of recent experiments. These two 
limits correspond respectively to paths of constant entropy or constant energy in the 
phase diagram; the final molecule fraction and condensate fraction depend strongly on 
which of these paths is taken. The interpretation of the recent experiments of Regal 
et al. Q] and Zwierlein et al. [3] requires specific consideration of the time dependence 
of the sweep. This is because the only temperature measurements reported are those 
taken before the onset of the sweep. These experiments appear to have been done 
in the adiabatic regime. Consequently it is crucial to plot theoretical predictions 
versus the initial temperature or entropy in an adiabatic phase diagram, rather than 
the more traditional one (e.g. Figure ^ used in previous interpretations of these 
experiments [3 Ej . Our simple model for a trapped system yields an adiabatic phase 
diagram in qualitative agreement with experimental data in the BEC limit, where a 
critical initial temperature of Tic ~ 0.2Tp is found. Here Tp is the Fermi temperature. 




Figure 2. Illustration of a Feshbach resonance. Panel A shows that a bound 
state in a closed channel can be in resonance with the scattering threshold of 
the open channel. Panel B illustrates a two-level model to describe the quantum 
dynamics of a pair of colliding atoms, where one level corresponds to a scattering 
state and the second level is the bound state of the closed channel. 



Before presenting our model and results, it is useful to place our work in the 
context of previous work done on this system. There are two main aspects of the 
problem we address that have been treated essentially as separate questions in the 
literature. The first question has to do with the formation of molecules as the 
resonance energy is varied in time [3 El El ^1 ^Q ^1 ^1 ^M ^1 and the second 
issue is the determination of the equilibrium phase diagram ^1 El HI HI El El HOI ■ 

A Feshbach resonance occurs between two atoms when the energy eres of a closed- 
channel bound state coincides with the relative energy of the atoms |23 1221 ESI- ^ 
diagram of a Feshbach resonance is shown in Figure |2K- The difference in magnetic 
moments between the open and closed channels allows the resonance to be tuned by 
adjusting an external magnetic field B. Due to the exchange interaction between 
valence electrons of the colliding pair of atoms, a coupling exists between the two spin 
channels. Due to this coupling, the pair of atoms spends a finite time in this closed 
channel bound state during the collision. This has a dramatic effect on the atomic 
interactions in the gas and provides a way for pairs of atoms to form a bound state, 
even in the absence of a third atom influencing the collision. 

Most of the work on the formation of Feshbach molecules is formulated to address 
highly nonequilibrium quantum dynamics that occur on a timescale much shorter than 
the time required for the system to relax to thermodynamic equilibrium. The starting 
point in such an approach is to first understand the process of molecule formation 
of an isolated pair of colliding atoms [3 |S] . One can treat the system as a two-level 
system, where one level corresponds to a superposition of scattering states in the open 
channel and the second level is the bound state of the closed channel, as illustrated 
in Figure EJB. The resulting coupled equations of motion can then be used to obtain 
a Landau-Zener adiabaticity criterion for establishing how slowly the _B-field must be 
ramped to achieve a full transfer to the bound state X i^ ■ More advanced treatments 
including many body effects have so far focussed on a gas of bosonic atoms near 
zero temperature ^■' 10l[mil2[ll3| . In these theories, the system is typically far from 
thermodynamic equilibrium. This treatment is justified, since most of the experiments 
on Feshbach resonances in Bose gases are forced to probe the system for a very short 



time due to inelastic collisions, which heat up the gas sample |24l 1251 1^ 1271 128| . 

In a Fermi gas, however, the same type of inelastic processes are suppressed by 
several orders of magnitude |^ OHl |^ , resulting in a very long lifetime of the gas 
sample |33 ESI 1^ IHH] ■ This permits experiments to be carried out over a time long 
enough for the system to relax to thermal and chemical equilibrium. In this case, a 
statistical description may be more appropriate than a microscopic description. The 
dynamics of molecule formation can then be treated using a kinetic theory approach 
with Boltzmann-type equations that describe the coUisional relaxation processes, 
which set the timescale for adiabaticity |36| . 

In the thermodynamic approach, a good starting point is to first map out the 
equilibrium phase diagram of a Fermi gas with resonant interactions, a qualitative 
sketch of which is shown in Figure 1. The first calculations of the phase diagram 
for this type of system were done by Eagles [23, Leggett [SH], and Nozieres and 
Schmitt-Rink 39' using a single channel model for a homogeneous Fermi gas. This 
work was improved upon later by Hausmann 00] and Chen et al. 0J. More 
recently, the phase diagram was calculated using a two-channel model for both 
homogeneous ^J El El El E] and trapped gases JHl 1211] ■ One of the main goals 
of all of these calculations is to determine the transition temperature separating the 
superfluid and normal phases, denoted by the blue line in Figure 1. The important 
question of how the system traverses from one point in the phase diagram to another 
as the resonance energy is varied was not addressed in these calculations. 

To address the question of molecule formation within equilibrium thermodynam- 
ics, one can consider the two limiting cases of a very slow sweep of Cros or a sudden 
jump. If ercs is varied at a rate much lower than the coUisional relaxation rate, the 
system will remain in equilibrium and the process will be adiabatic |36j . In this case 
the entropy is conserved during the sweep. Some aspects of this case were explored 
by Kokkelmans et al. for a normal gas W? and by Carr et al. in the superfluid 
regime .43, ,^]. In the other limit of changing eros suddenly, on a timescale much 
shorter than the relaxation time, the entropy increases. In this case, the energy is 
conserved as the system relaxes to a new temperature. Chin and Grimm discussed 
this case at length, treating an ideal gas mixture of atoms and molecules in the classical 
regime T>rc gl]. 

In this paper we present a detailed study of both cases of adiabatic and sudden 
variations in Cres- Our work extends the work of Chin and Grimm into the quantum 
degenerate regime. Our simple model, which neglects interactions, describes BEC of 
molecules but does not allow for BCS pairing superfluidity on the eres > side of the 
resonance. The primary goal of our work is to develop a conceptual understanding 
of molecule formation in the adiabatic and sudden limits within the framework of 
statistical mechanics. For this purpose the ideal gas is an excellent toy model. 

Our paper is outlined as follows: in Section II we present the equilibrium theory 
for an ideal gas mixture of fermionic atoms and bosonic molecules that are in thermal 
and chemical equilibrium. Using this formulation, we calculate the phase diagrams 
in Section III, showing the molecule fraction and the fraction of condensed molecules 
versus the temperature T and the resonant energy eies- In IIIA we consider the 
adiabatic limit, and in IIIB the sudden case. We then close with some concluding 
remarks in section IV. 



2. Equilibrium theory 

We consider an ideal gas mixture of (pseudo) spin- 1/2 fermionic atoms and zero spin 
bosonic molecules in an anisotropic harmonic trap with frequencies tOx, tOy, and luz 
along the respective axes. We assume the trapping force on the molecules is twice 
that of the force on the atoms, so that the atomic and molecular frequencies are 
identical gUJ . The equilibrium distributions for the atoms and molecules are given by 
Fermi-Dirac and Bose-Einstein distributions, respectively 

^''^'^' " e(e-M.)/fcBT„ + 1' ^^' 

•^™(^) " e(<!+e„=-/^™)/'cBT,„ _i- (2) 

where a G {|, J,} signifies the spin state of the atoms and fee is Boltzmann's constant. 
In this work we restrict our attention to an equal spin mixture /i| = /ij^ = /ia and 
T| = Tj = Ta, so that /| = /x = fa- We assume the atoms and molecules are in 
thermal equilibrium so that the temperatures are equal Ta = Tm = T. We also assume 
that the system has relaxed to chemical equilibrium due to collisions that convert pairs 
of atoms into molecules and vice versa |3()l I45| . In this case /im = 2/^^ = 2/i. The 
energy of a molecule is displaced relative to that of the atoms by eros, the magnitude 
and sign of which can be adjusted. 

To obtain the phase diagrams, we need to calculate the population, energy, 
and entropy for each component. The atom population Na{T,fi) (in a single spin 
component) and noncondensed molecule population Nm(T,fi, ercs) are 



pW = o7^-W' (5) 



NaiT,fi) = / dep(e)/a(e), (3) 

Jo 

/>oo 

N^{T,ti,e,,s)= / dep{e)U{e)- (4) 

Jo 

The density of states in a harmonic trap is 

where w = {ojxtiJyUJz)^'^- Our semi-classical formulation, which replaces the sum over 
discrete states by an integral, is strictly valid when k^T ^ huj. When this inequality 
is not satisfied, shell effects in the Fermi sea can play a role 07|SH1- The integrals can 
be evaluated and expressed in terms of Fermi ^„(z) and Bose Qn{z) integrals defined 

by 

where T{n) is the gamma function. The Bose function is equal to the Riemann-zeta 
function C(") when z = 1, i.e. Gn{z = 1) = C{n). 

In terms of these functions, the populations are 

Na{T,^i) =(^)V3(Za), (8) 

Nra{T,fi,e,,s) = {^yOsiZn,), (9) 



where the atom and molecule fugacities are defined by 

Za{T,fJ.) 



= gP/^BT 



The atom energy Ea{T,ii) and noncondensed molecule energy E„i{T, ii,ercs) are 



(10) 
(11) 
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= SNakBT 



^A{Za) 
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= iNmkBT 



^■i{Za) 

dep(e)(e + eros)/m(e) 

Qi{Zn 



(12) 



+ e^csNr, 
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Finally, the atom entropy Sa{T, /i) and noncondensed molecule entropy Sm{T, /i, e 



(13) 



are 



/•OO 

Sa = -k-B, I dep{t) [fa In /„ 
Jo 

.Pi{Za) M 



+ (l-/a)ln(l-/,) 



k^Na 



knT 



H{Za) 

Sm= - ks dep{e) [/„ In /„ - (1 + /,„) ln(l + /„)] 

Jo 

C4(z„) (2^-ercs)" 



= /cB^m 



(14) 



(15) 



G?,{Zm) k^T 

At a low enough temperature, the molecules form a Bose-Einstein condensate. 
Above the transition temperature T > Tc(ercs), fJ- < ercs/2 and the population of 
condensed molecules is zero. In this case the total population of atoms is 

A^tot(T, fi, e,es) = 2A^,(T, ^i) + 2iV„.(T, fi, e.-e.), (16) 

where the factor of 2 in the first term describes the two spin components and in the 
second term there are two atoms in each molecule. For a given value of the population 
A'^, the chemical potential is obtained at a specified temperature T and resonance 
energy e^es by solving for /i in the equation 

^tot(T,Ai,e,es) = ^. (17) 

Below Tc(eres), fJ- — ercs/2, so that z„i = 1. The condensate population Nmc{T, Cios) is 

N„UT, ercs) = N/2 - Ntot{T, fi = e,,,/2, e,,,)/2. (18) 

Below Tc the noncondensed fraction N^n reduces to the expression 



A"„i(T, eres/2,er. 



«^)(^)- 



(19) 



The energy of condensed molecules is ercsNmc{T, Crcs) and the condensate entropy is 
zero. 

On the positive detuning side of the resonance eics > 0, we can obtain an analytic 
expression for Tc by approximating the population of atoms by the T = limit |49| 



A^a(r = 0,M = ercs/2) 



6\2hLu/ 



(20) 
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Substituting Eq. (fT^ and Eq. l(^ into Eq. (^Hj), setting N^c — 0, and solving for 
T = Tr, we obtain 






^tot _ 1 /£rcs Y 

2 6V2nw/ . 



1/3 



(21) 



C(3)i/3fc, 

Taking Tc = in Eq. (|21|l . we can solve for the critical value Ec, below which 
condensation occurs at zero temperature 

e, = 2fcBTF. (22) 

The Fermi temperature Tp in a trap is defined by /cbTf = {SNY^^huj. On the negative 
detuning side ercs < 0, in the limit lerosl/fceT^ ^ 1, we may set Na = in Eq. 118|) . 
Solving for Tc we obtain the expected result ^l§^ 

^ = [6C(3)]-i/3« 0.518. (23) 

J-F 

Equations (|22|l and (|23|l define the right and top boundaries of the condensation 
region in the phase diagram for an ideal gas mixture of fermionic atoms and bosonic 
molecules. 

3. Phase diagrams 

We now turn to the numerical calculation of the molecule fraction ?7m(eros,T) and 
fraction of condensed molecules ?ymc (crcs, T) defined by 

^™ = ^, (24) 

Vmc — -jr; J (25) 

where Nm — N„i + Nmc is the total molecule population. In our calculations, we fix 
the total population N . Then, for a given value of T and eios, we numerically solve 
for the value of /i that satisfies Eq. (|16|l . As we approach the condensation region, 
the chemical potential /i — > eros/2. Once we cross into the condensation region, the 
condensate population is calculated according to Eq. lfTF|l . 

The results of the calculations are shown in Figure O We plot the molecule 
fraction r]m{(-rcs,T) in Figure |3J\ versus T and ercs, where dark blue corresponds to 
no molecules in the gas and dark red means there are only molecules. We also plot 
the condensate fraction rjmcisics, T) in Figurc|33. The energy scale is set by the Fermi 
energy k^Tp. By normalizing both axes by this energy, the resulting phase diagrams 
for an ideal gas mixture of harmonically trapped atoms and molecules are " universal" , 
in the sense that they look identical for any set of physical parameters used in the 
calculation (i.e. population N , mass m, and trap frequency uj) . Our calculation 
of the condensate fraction agrees with the prediction that the condensation region is 
bounded on the right by Cc = ^k-oT-p and on the top by Tc ~ 0.518Tf. 

In Figure 13'^, the green region signifies a dissociation region where half of the 
atoms have been converted to molecules. We define edis(r) as the value of eics where 
rjm = 0.5. In recent experiments, this dissociation region is associated with the position 
of the resonance 32, 33 . Our calculations for an ideal gas demonstrate that e^isiT) 
depends on the temperature. We can obtain an analytic expression for edis(T') in 
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Figure 3. Phase diagrams. The molecule fraction »?m(eres,T) is plotted in panel 
A and the condensate fraction »7mc(eres,T) is shown in panel B. The color scale 
is shown on the right. In panel A, the high (solid line) and low (dashed line) 
temperature limits of edis(T') are shown. 



the two limits of low or high temperature. In the low temperature limit, we use the 
Sommcrfeld expansion of Tn{e^) for x » 1 [501 
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Using Eq. (PS|l . we can solve for Cdis in the equation Na{T, edis/2) = A^/4 

24/3^2 y2 



Cdi 



,{T « rF)/fcB = 22A''rF - 



(26) 



(27) 



3 Tf 

In the limit of T ^ Tf , the atom and molecule distributions go over to the Maxwell- 
Boltzmann distributions and it is straightforward to verify that 

edis(T » rF)/fcB = - [3 ln(r/rF) + In 12] T. (28) 

The high and low temperature limits of edis(T') are shown in Figure|31by the solid and 
dashed black lines, respectively. 

To address the question of which path is taken in the phase diagram as the 
resonance energy is varied, we consider the two limiting cases of an adiabatic sweep 
and a sudden change in Cros- 



3.1. Adiabatic sweep of e^cs 

If the resonance energy is varied on a timescale much longer than the time needed 
for the gas to coUisionally relax to equilibrium, then the system will follow a path of 
constant entropy in the phase diagram. The total entropy of the gas is given by 

'S'tot {T, H, eres) = '2Sa (T, fj.) + S„i (T, /X, Eres) ■ (29) 




-4 -2 

resonance energy 



4 



4-2024 
resonance energy f^Tcs/^BTY 



Figure 4. Contours of constant entropy. Each solid white line defines the path 
followed by an adiabatic sweep of the resonance energy eres- For clarity, we plot 
the contours on top of the phase diagrams with the same color scale as in Figure 
3. The dashed line corresponds to an initial temperature of Ti/T-p = 0.1. 



Note that the condensed molecules do not contribute to the entropy in the system. 
Starting on the right side of the phase diagrams at some initial temperature T^ and a 
large positive value of eros: there are no molecules and the initial entropy is given by 



S tot {Tt,fii) = kBN 



fJ-i 



(30) 



The resonance energy Crcs is slowly lowered to a new value Crcs,/- The new equilibrium 
solution is obtained by solving for Tf and /i/ that satisfy the constraints of the total 
population and entropy being conserved. Above Tc this corresponds to solving for Tf 
and fj,f in the equations 

Ntot{Tf,flf,£rcsj)^N, (31) 

Stot{Tf,flf,£rcsj) = Stot{Ti,fJ.i). (32) 

Below Tc, fif = ercs,//2, so only Eg. IS^ must be solved. 

In Figure ^ we plot the contours of constant entropy (white lines) overlaying 
the phase diagrams. Each line defines the path followed by an adiabatic sweep of the 
resonance energy Crcs- Starting on the right side at a large value of eics, the temperature 
increases as eics is lowered adiabatically, but eventually reaches a constant value when 
all of the atom pairs have been converted into molecules. We can understand the 
temperature increase intuitively: as pairs of atoms are converted into molecules, the 
system loses degrees of freedom, so the molecules must heat up to conserve entropy. 
Because an adiabatic process is reversible, it is possible to cool the gas by starting on 
the left side with all molecules and increasing ercs- This useful idea was first pointed 
out by Carr et al. as a means to achieve a low enough temperature to get into the 
BCS region 03). We emphasize that all of the atoms can be transfered into molecules 
if the sweep is done adiabatically. 

An alternative way to view the same information in Figure^is to plot rjm and rj^nc 
versus eres and the entropy S, which we show in Figure [S] We label the y-axes by the 
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Figure 5. Adiabatic phase diagrams. The molecule fraction r;m(eros, S) is plotted 
in panel A and the condensate fraction T;mc(eres, S) is shown in panel B. We label 
the jy-axis, which corresponds to the entropy, by the initial temperature starting 
out with all atoms. The dashed white line corresponds to the dashed line of 
constant entropy in Figure |1] 



initial temperature Tj starting with all atoms; the corresponding value of the initial 
entropy can be obtained from Eq. (|30|l (where fii is chosen to give a total population 
of N). The resulting adiabatic phase diagrams look qualitatively different from the 
phase diagrams shown in Figure|3 In recent experiments, an adiabatic phase diagram 
was obtained by measuring the condensed molecule fraction after the magnetic field 
was adiabatically lowered starting with all atoms at some initial temperature T^ mS] . 
It would be interesting to repeat these experiments and measure the final temperature 
of the gas after lowering the magnetic field, and thus obtain a phase diagram like the 
one shown in Figure 13 

The dashed line in Figure |31 corresponds to the dashed contour shown in Figure 21 
for an initial temperature of Ti/Tp — 0.1. In Figure we plot the molecule fraction 
and condensate fraction along this constant entropy contour. As eres is lowered, the 
condensate fraction rises sharply and then decreases to some final value rjmcj after 
all the atoms have been converted to molecules. This same qualitative behavior was 
observed in recent experiments 0E|- 

We can obtain analytic expressions in the high and low temperature limits for 
the final temperature and condensate fraction after all the atoms are converted to 
molecules. The initial entropy is Stot{Ti) = 2Sa{Ti, fii), where Sa is given in Eq. ifT^ . 
Using the Sommerfeld expansion in the low temperature limit, we find 

H 



Stot{T^)^TT^kBN^. 



(33) 



The final entropy is Stot(Tf) 
taking fif = eros/2, we find 



~ S'„i(T/, ei.cs/2,ercs)- Using expression Eq. ((T3|) and 
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t(T/)-4C(4)fcB( 
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(34) 
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Figure 6. Molecule fraction J7m(T, tros) (solid red line) and condensate fraction 
Vmc{T, ercs) (blue dashed line) for an adiabatic ramp. These curves correspond to 
the dashed lines in Figures f^ and 151 where the initial temperature is T/Tp = 0.1 
starting with all atoms. 



Solving for Tf in the equation StotiTf) = S'tot(T'i) and using N — (k-BTF/huj)^ /'S, the 
final temperature in the low temperature limit T/Tp <C 1 is 



Tf ^ / TT^ \ 1/3 / T, \ 1/3 

Tf V12C(4)/ \Tf. 

nl/3 



(35) 



The result Tf oc T^'"^ agrees with the results of Carr et al. ^^ . Their model treats 
mean field interactions in the condensate, so the constant of proportionality is different 
in the two models. 

In the high temperature linut T/Tp 3> 1, .7-n(z) w z and Gn{z) w z. Using this in 
the expressions for the entropy Eq. (|14|) and Eq. (|15|l . and using Eq. ||SJ| and Eq. 
to eliminate the chemical potentials, we obtain the initial and final entropies in the 
high temperature limit 

hui \N-\ 



S tot [Ti) 
StotiTf) 



kBNH 

1. 



I LV/crT,/ 2 Ji' 



-knNU-\n 



ksTj 2 
KkBTfJY 



(36) 
(37) 



Solving StotiTf) 
T/Tf » 1 is 



StotiTi), the final temperature in the high temperature limit 



^=6^/3^^/^''— V 



In FigureElwe plot the final temperature as a function of the initial temperature. 
The black dotted line was obtained from a numerical solution. The solid red line 
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Figure 7. Final temperature versus tlie initial temperature after all the atoms 
have been converted to molecules. 
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Figure 8. Condensate fraction rjmc.f after all the atoms have been converted to 
molecules versus the initial temperature. 



corresponds to the high temperature hmit Eq. (|38|l and the blue dashed Hne to the 
low temperature limit Eq. (|35|l . Both approximations agree well with the numerical 
solution in their appropriate regions of vahdity. 
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We can calculate the condensate fraction at the end of the sweep in the hmit of 
l^Tcsjl/kBTp ^ 1 as a function of the intial temperature using the low temperature 
result in Eq. H35|) . In terms of the final temperature, the condensate fraction is given 
by r]rnc,f = 1 - (Tf/Tcf, where T^ = [6C{3)]-'^/^Tf. Substituting T/ from Eq. ^. 
the condensate fraction in this limit varies linearly with the initial temperature as 

V ,-i_!l!MZl (39) 

The critical initial temperature predicted by this equation, above which no condensed 
molecules exist when all the atoms have been converted to molecules, is Tj^c/Tp « 
0.183. In Figure IHl we plot the condensate fraction after an adiabatic sweep versus 
the initial temperature. We compare the low T approximation to the full numerical 
solution. The deviation near Ti^c arises from using only the first two terms in the 
Sommerfeld expansion in the low T approximation. 

Our simple model calculation helps to clarify the interpretation of experimental 
data concerning the critical temperature for superfluidity. In the BEG limit of negative 
detuning, the critical temperature observed in recent experiments is approximately 
0.2Tf Q 01' At first sight, this seems in contradiction with the prediction of 
Tc/Tp = 0.518 for a trapped Bose gas ^^1 and seems more consistent with the 
homogeneous gas result [SJ |H] of Tc/Tp = 0.218. However, when one realizes that 
the experimental data is plotted versus the initial temperature, the observed value of 
Tc/Tp « 0.2 is consistent with the expected result for the critical initial temperature 
Ti^c for a trapped ideal gas with an adiabatic sweep, as shown in Figure |H1 and in the 
left side of Figure O?. 

3. 2. Sudden change of eres 

In the previous section we considered the situation where ercs is varied on a timescale 
much longer than the coUisional relaxation time for the gas to equilibrate. In this 
adiabatic limit, the gas dynamically follows along the constant entropy contours shown 
in Figure 01 and the process is reversible. We now consider the opposite limit of a 
sudden change in Cres- In this case, it is the energy that is conserved rather than the 
entropy, and the process is not reversible. In the two sections below, we consider two 
different initial states, starting out with either a gas of only atoms (the right side of 
the phase diagram) or with a gas of only molecules (the left side of the phase diagram). 

3.2.1. Starting with all atoms Suppose the gas of fermionic atoms is cooled down 
to some temperature Ti and that ercs.i/^B^p 3> 1, so that there are no molecules 
initially, (/ii must be chosen to satisfy Na{Ti,iJ,i) — N/2). The resonance energy is 
then suddenly lowered to a new value ei-csj over a time much shorter than all coUisional 
relaxation timescales. The initial total energy of the gas immediately following this 
jump in Ercs is 

Etot{T„fi,) = 2Ea{T,,fii), (40) 

where Ea{T,iJ,) is defined in Eq. (|13(l . The system subsequently relaxes from this 
initial state to a new equilibrium state with a different temperature Tf and chemical 
potential /x^, during which time pairs of atoms may form molecules. The final energy 
is 

EtotiTf,flf,ercsj) = 2Ea{Tf,^J,f) + Ern{Tf,^J,f,ei-cs,f) 

+ eres,fNmc{Tf,eresj)- (41) 
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resonance energy fres/^BTp resonance energy ii-cs/kBTF 

Figure 9. Lines of constant energy starting with all atoms initially. The initial 
temperature for the lines start at T^ = and increase in increments of O.lTp. 
As in Figure |1] we plot the contours on top of the phase diagrams showing the 
molecule fraction in panel A and the fraction of condensed molecules in panel B. 
For clarity, we show hypothetical intial and final points in panel A. 



Because the Haniiltonian for the system is conservative and does not depend on time 
after the sudden jump in e,cs, the final energy must be equal to the initial energy 

£'tot(T'/,M/,eros,/) ^ Etot{Ti,^j,i). (42) 

The values of Tj and /i/ are chosen to satisfy Eq. 142|) and number conservation. 

In Figure El we plot the lines of constant energy starting with all atoms initially. 
The initial temperature of the lines (on the right side of the plots) increases in steps of 
O.lTp, starting with the lowest line at zero temperature. The region below the lowest 
line is a forbidden region of the phase diagram that can not be reached by suddenly 
lowering eics- In Figure EK we show a hypothetical final point starting at an intial 
temperature of T^ = 0.4Tf. We emphasize that the system does not dynamically follow 
along the path between the initial and final points, in contrast to the adiabatic case. 
After the sudden decrease in eresi the system evolves out of equilibrium away from the 
initial state and then relaxes to the final equilibrium state at a new temperature Tf. 

From Figure EK we can see that it is impossible to transfer all of the atoms to 
molecules by suddenly decreasing eics- This property was also noted by Chin and 
Grimm in Ref. ;45[ for a classical gas. As ej-es decreases, the molecule fraction rjjn a-nd 
condensate fraction rjmc each reach a maximum and then decrease. These maximum 
values max{?7„i} and max{77mc} increase as the initial temperature Ti is lowered. In 
Figure E| we plot ijm (solid red line) and ijmc (blue dashed line) versus eres for an 
initial temperature of Ti — 0. At a value of Cros = 2kBTF, all of the molecules are 
condensed max{77mc} = 1- As ercs is lowered, the temperature T increases and the 
condensed fraction decreases sharply to zero. The molecule fraction increases to a 
maximum value of maxlr^m} ~ 0.7, then decreases slowly. 



3.2.2. Starting with all molecules We now consider the situation where we start out 
on the negative side of eres with a gas of bosonic molecules at some initial temperature 
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Figure 10. Molecule fraction r)m{T, tros) (solid red line) and condensate fraction 
rimc{T,tics) (blue dashed line) for a sudden ramp. These lines correspond to the 
bottom line shown in Figure |5] which starts at an initial temperature of Ti = 0. 



Ti and that lercs.il/fcBTp ^ 1, so that there are no atoms initially, {^i must be chosen 
to satisfy Nrn(Ti,^j,i,ercs,i) + Nmc{Tt, fj,t,eres,i) = N/2). The resonance energy is then 
suddenly raised to a new value Cres,/ over a time much shorter than all coUisional 
relaxation timescales. The initial total energy of the gas immediately following this 
jump in ercs is 



EtotyTi, ^ii, Cros,?:) — '^Nm{Ti, fli 



)knTi 



Gi{Zr, 



GsiZn 



N 

+ ^ErcsJ, (43) 



where Zm,i = e^'+''^'==' is the initial molecule fugacity. The first term describes the 
initial thermal energy of the gas. Note that it is the final value of eics./ appearing in the 
second term describing the internal molecular energy. The system subsequently relaxes 
to a new equilibrium state with a different temperature Tf and chemical potential /jj , 
during which time some of the molecules may dissociate into pairs of atoms. The 
final energy is given by Eq. H41() . The final values of Tf and /// are chosen to satisfy 
Etot{Tf,^f,e,-csj) = Etot{Ti,^i,e^cs,i) and number conservation. 

In Figure ^J we plot the lines of constant energy starting with all molecules 
initially. The initial temperature of the lines increases in steps of 0.2rF, starting with 
the lowest line at zero temperature. The region below the lowest line is a forbidden 
region of the phase diagram that can not be reached by suddenly raising e-^cs- In the 
previous case of starting with all atoms, there is an energy barrier preventing a full 
transfer of all the atoms into molecules. Here we find that all of the molecules can 
dissociate as Errs is increased. 
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Figure 11. Lines of constant energy starting with all molecules initially. The 
initial temperature for the lines start at T; = and increase in increments of 
0.2Tp. The starting points of the upper lines lie out of the range shown in the 
graph. 



4. Conclusions 



In this paper, we have studied the thermodynamics of molecule formation in an ideal 
gas mixture of fermionic atoms and bosonic molecules that are in chemical and thermal 
equilibrium, where the internal energy of the molecules can be adjusted relative to that 
of the atoms. We considered both limiting cases of sweeping the resonance energy 
adiabatically or changing it suddenly. Our work is the natural extension into the 
quantum degenerate regime of the earlier work by Chin and Grimm ^45, who studied 
a classical gas and has some overlap with the recent work of Carr et al. |43[ E| . 

Our toy model serves as a useful guide for studying the problem of molecule 
formation in the adiabatic and sudden limits. In particular, we have given a 
prescription for calculating the adiabatic phase diagrams, shown in Figure [S] in 
which the molecule fraction and fraction of condensed molecules are plotted versus 
the resonance position Ej-cs and entropy S. These should be contrasted to the phase 
diagrams shown in Figurc|3 where quantities are plotted versus temperature T rather 
than entropy. In the experiments of Regal et al. T and Zwierlein et al. [H], the 
condensate fraction was measured as a function of the initial temperature and the 
magnetic field controling the resonance. The resulting graphs in Figure 4 of Ref. PP 
and Figure 5 of Ref. T are adiabatic phase diagrams, analogous to Figure |SJ3; they 
are not the more intuitive phase diagram shown in Figure |2|3 jSJ |H] . This distinction 
is important to realize when comparing theory to experiment. 

As an example, the theory of a trapped gas previously developed by Ohashi 
and Griffin jJS] predicted a BEG transition temperature of T^ ~ O.STp. Subsequent 
treatments of a homogeneous gas by Falco and Stoof |S] and Diener and Ho [H] obtained 
a transition temperature of Tc ~ Q.2T-p. However, none of these theoretical values 
can be compared directly with any temperature that was actually measured in the 
experiments of Regal et al. 1^ and Zwierlein et al because the only temperatures 
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reported there are those of the Fermi gas at the start of the magnetic field sweep. 
The construction of the adiabatic phase diagram that is appropriate to compare our 
model with experimental data shows that an atomic gas with initial temperature 
of Ti « 0.2Tf is adiabatically converted to a molecular gas with a temperature of 
Tf w 0.5Tp, which is equal to the critical temperature for a trapped ideal Bose gas. 

In order to make a detailed comparison with experiments, a more advanced 
many-body theory that properly treats the effect of the resonant interactions must be 
used in order to describe the crossover to BCS superfluidity in the positive detuning 
region 1161 1171 1511^ 1181 [TO| . Furthermore, in experimental data the resonance position 
is typically expressed in terms of the external magnetic field B rather than the energy 
of the resonant state eros(^); this must also be taken into account when comparing 
theory to experiment. 
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